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Myelin figures are interfacial structures formed when certain surfactants swell in excess water. 
Here, I present data and model calculations suggesting myelin formation and growth is due to the 
fluid flow of surfactant, driven by the hydration gradient at the dry surfactant /water interface; a 
simple model based on this idea qualitatively reproduces the various myelin growth behaviors ob- 
served in different experiments. From a detailed experimental observation of how myelins develop 
from a planar precursor structure, I identify a mechanical instability that may underlie myelin for- 
mation. These results indicate the mixed mechanical character of the surfactant lamellar structure, 
where fluid and elastic properties coexist, is what enables the formation and growth of myelins. 



PACS numbers: 82.70.Uv, 61. 30. St, 62.20.D- 

I. INTRODUCTION 

A bit of soap left in a wet dish turns into a gelatinous 
mass — the sight may be unexceptional, but the micro- 
scopic structure of wet soap (hydrated surfactant) can 
be exceptionally diverse. The disparate interactions of 
the surfactant molecules with themselves and with wa- 
ter leads to the formation of aggregate structures, from 
spheres to cylinders (spherical and worm-like micelles) to 
extended membranes (bilayers); these can further orga- 
nize into yet more complex structured phases [1]. Here, 
temperature, pressure, and concentrations — of surfac- 
tants, solvents, salts, etc. — are the principle coordinates 
of a potentially very complex phase diagram. 

Research into lyotropic (concentration-dependent) sur- 
factant/water phases have largely focused on the equi- 
librium phase diagram. The kinetic processes through 
which the non-equilibrium initial state of pure surfactant 
and pure water arrives at equilibrium has been subjected 
to relatively few studies. Due to the strong concentration 
gradients, the interface between pure surfactant and pure 
water can consist of multiple domains of different phases, 
as well as metastable structures whose origins are not un- 
derstood, evolving in spatially and temporally complex 
ways [21 [11 m |5] . But even in surfactants whose phase di- 
agram is fairly simple, the non-equilibrium kinetics can 
result in striking interfacial instabilities such as myelin 
figures. 

First described in 1854 by Rudolf Virchow, who pio- 
neered the use of microscopy in clinical pathology [3], 
myelins are /im-diameter cylinders consisting of con- 
centrically stacked bilayers, and are formed when some 
poorly-soluble surfactants swell in excess water [51 [31 [7J [HI 
[HI Uni HH [H] • Their sometimes dramatic growth (reach- 
ing lengths ~ 100 /xm in a few seconds), together with 
their sinuous figure, give myelins a strangely life-like ap- 
pearance. In the 150 years since their discovery, the 
origin of myelins have remained mysterious. The equi- 
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FIG. 1; Myelin flgures viewed using dark field microscopy. 



librium structure of myelin- forming surfactants consists 
of planar surfactant bilayers stacked one on top of an- 
other, separated by thin layers of water. Myelins, with 
tightly curled bilayers, have large bending elastic ener- 
gies. Why then, do these seemingly energetically unfa- 
vorable structures form? Superficially, myelin formation, 
with the growth of finger-like structures from an initially 
flat interface, is reminiscent of non-equilibrium fingering 
processes such as the MuUins-Sekerka instability of a so- 
lidification front [13]. But a successful model for myelin 
formation based on such instabilities has yet to be intro- 
duced. 

The present paper proposes an alternative scenario. 
Due to the large concentration gradients at the concen- 
trated surfactant/ water interface, strong currents (of wa- 
ter and/or surfactant) must be present; these can focus 
large stresses at the interface, possibly leading to me- 
chanical instabilities. I suggest myelin formation is due to 
such a mechanical instability. Interestingly, because the 
mechanical properties of the surfactant lamellar structure 
is anisotropic — depending on the direction, it behaves 
either as a fluid or as an elastic solid, the myelin instabil- 
ity is one that has both fluid and elastic characteristics. 

This paper is organized as follows. Section II describes 
three distinct myelin growth experiments; their results 
are complementary, and together brings a fuller, more 
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FIG. 2: (color online) The contact experiment. Surfactant is 
sandwiched between glass slides, and water (from a syringe) 
is introduced into the gap via capillary suction. Myelins form 
rapidly along the interface upon contact between water and 
surfactant. 

detailed picture of how myelins form and grow. In Sec- 
tion III, motivated by the varying myelin growth behav- 
iors observed in these experiments, I introduce a simple 
model for myelin growth based on the fluid flow of surfac- 
tant. The driving force behind this flow is the hydration 
gradient at the surfactant /water interface. Compared 
with experiment, this model qualitatively reproduces the 
observed myelin growth behaviors. Finally, Section IV 
describes in detail the formation of a myelin from a pla- 
nar lamellar structure as observed in one of the experi- 
ments. I identity a mechanical instability which may un- 
derly myelin formation, and which gives good agreement 
with experimentally observed morphology and dynam- 
ics. Here the coexistence of fluid and elastic properties 
in the surfactant lamellar structure is crucial: in order 
for a myelin to form, the surfactant lamellar structure 
has to flow as a fluid and buckle as an elastic solid. 



II. EXPERIMENTAL METHODS 

Some of the results from the experiments described be- 
low have been previously reported |12j : briefly summa- 
rized, these are: (1) myelins form laterally from the edges 
of their parent lamellar structures; (2) myelin growth is 
due to the flow of surfactant from the base and follow 
different time-dependences in different experiments; and 
(3) myelins grow and persist only in the presence of ther- 
modynamic or mechanical stress. In Sections III and 
IV, I will revisit these results, present additional data, 
and compare experimental results with theoretical mod- 
els. For this Section, I will restrict myself to describing 
the details of various experimental methods. 




FIG. 3: (color online) (a) The immersion-and-puncture ex- 
periment as observed under the microscope, (b) Schematic 
illustration of the immersion-and-puncture experiment; pre- 
sumably, the puncture as a hole of some effective size ro which 
eases water invasion into the surfactant plaque. 



A. Contact Experiment 

The simplest method to produce myelins to bring into 
contact a mass of appropriate surfactant with water. The 
surfactant I use is dilauroyl phospaditylcholine (DLPC), 
a common phospholipid. DLPC is dissolved in chloro- 
form (concentration: 20 mg/mL); a drop of the solution 
is cast onto a clean glass surface, and dried in dry nitro- 
gen overnight. Just before use, the surfactant plaque is 
softened by brief exposure to the humidity from a warm 
water bath. The softened surfactant is carefully scrapped 
off and sandwiched between two glass slides, using a thin 
(0.15 mm) glass coverslip as spacer. Water is then in- 
troduced into the gap via capillary action, and contacts 
the surfactant plaque along an extended interface. Upon 
contact, myelins in dense bundles form and rapidly grow 
along the water/surfactant interface (Fig. |2]). 

To measure myelin growth, the experiment is observed 
via digital video-equipped dark-field microscopy in low 
magnification. The long thin myelins strongly scatter 
the indirect illumination, and appears bright while the 
bulk mass of surfactant remains dark. Thus the bright- 
ness of the video field is proportional to the amount of 
myelins formed. From the the average brightness of the 
myelin bundles, the average length of the myelins can be 
computed. This allows myelin growth to be monitored 
with a much higher temporal resolution than is practical 
by manual measurement. 



B. Immersion-and-Puncture Experiment 

Contact experiments as described previously do not 
in general have a good control over the structure of the 
concentrated surfactant plaque. Ideally, one would like 
to start with a well-ordered lamellar initial structure and 
observe if and how myelins develop. 

To prepared the ordered initial structure, DLPC in 
chloroform is cast onto a clean glass-bottomed petri dish. 
Once the solvent has evaporated, the resultant surfactant 
plaques are incubated at 60° C and 100% relative humid- 
ity (r.h.) for up to two days. At the end of this process. 
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the plaques anneal to form well-ordered planar structures 
with mm-size domains. The plaques are cooled and dehy- 
drated to room temperature at 30% r.h. The dish holding 
the surfactant plaques is then carefully flooded with wa- 
ter. Surprisingly, on immersion, few myelins are formed 
from the well-ordered surfactant plaques; this is in con- 
trast to the dense thickets of myelins form on the surfac- 
tant/water interface in contact experiments. 

However, myelins can be induced to form by punc- 
turing the planar top surface of the surfactant plaques 
with a sharp needle (I use a tungsten or stainless steel 
electrophysiology electrode mounted on a micromanipu- 
lator). On pressing the needle into the plaque, one or 
more myelins will usually form at the site of the punc- 
ture and slowly grow in length. Presumably the punc- 
ture facilitates water invasion into the surfactant plaque 
(Fig. [3]). As a clear illustration of the lateral fluidity 
of the surfactant lamellar structure, by moving the nee- 
dle laterally across the plaque, the myelin can be readily 
dragged along without structural disruption. For earlier 
times, the growing myelin "crawls" along the surface, and 
its contour can be easily traced by a spline curve. As the 
myelin grows longer, its contour becomes convoluted and 
no longer quasi- 2D; at this point the myelin's length can 
no longer be reliably measured. 



C. Drying Drop Experiment 

The drying drop experiment is another way to grow 
myelins in a controlled way from a well-ordered planar 
initial structure. The surfactant used here is dimyristol 
phospaditylcholine (DMPC); it is identical to DLPC ex- 
cept it has alkane tails groups that are two carbon atoms 
longer. Dry DMPC powder is dissolved in ethanol at a 
concentration of 50 mg/mL. A small quantity (10 fiL) 
of the DMPC solution is pipetted into a microcentrifuge 
tube, then rapidly and throughly mixed with deionized 
water to form a vesicle suspension with final concentra- 
tion of 0.5 mg DMPC/mL. This method consistently pro- 
duces a fine suspension with a mean vesicle diameter of 

50 nm as determined by dynamic light scattering. The 
small amount of ethanol in the suspension does not ap- 
pear to influence myelin formation in these experiments. 
The same phenomenon can be observed in experiments 
using ethanol-free vesicle suspensions produced (labori- 
ously) by multiple extrusions across porous membranes. 

The experimental cell is constructed by epoxying a 
teflon o-ring on to a glass slide (Fig. |4|. A small drop 
of the DMPC suspension is placed inside the cell, and a 
glass cover slip is placed on top of the o-ring. The con- 
tact between the cover slip and the o-ring is intentionally 
left unsealed, allowing the drop to evaporate slowly: a 20 
IjL drop will take about 2 hours to dry completely. This 
geometry is best suited for observation using an inverted 
microscope; it can also be inverted, with the drop being 
pendant rather than recumbent, for observation using an 
upright microscope. The experiment is kept at 30° C. 




glass slide 



FIG. 4: (color online) Drying drop experiment: myelins grow 
from well-ordered lamellar pancakes located at the contact 
line of a slowly evaporating water drop. 



The experiment works via the "Coffee Stain Effect" 
dH US]. As the drop dries, the suspended surfactant is 
deposited by capillary flow at the pinned contact line as 
a ring-stain. The stain is mostly disordered, but small, 
well-ordered multilamellar "pancakes" are also formed 
right on the contact line. These pancakes grow slowly 
in size, and if one grows sufficiently large to invade back 
into the drop, it will form one or more myelins. 

Drop evaporation can be halted midway through an 
experiment, by injecting a small drop of silicon oil near 
the o-ring/cover slip contact using a flne pipette. Cap- 
illary action will pull the oil into the o-ring/cover slip 
contact, forming an air-tight seal. Once the sealed cell 
equilibrates, in ^ 1 min, myelins stop growing and resorb 
into the pancakes from which they had developed earlier. 

For this experiment to work effectively, it is important 
that the contact-line of the drying drop is pinned. The 
surface characteristics of commercially available glass 
slides can vary considerbly, but the following treatment is 
found to consistently yield the desired properties. Glass 
slides are sonicated sequentially in toluene, methanol, 
and deionized water. They are then etched in a reactive 
ion etcher using CF4, removing 50 nm of material, fol- 
lowed by a short etch with O2. Finally they are "aged" 
in dry nitrogen for 10-12 hours before use. 

Experiments are observed using transmitted bright 
fleld microscopy. Real-time videos and still images are 
captured directly to computer using a CCD camera at- 
tached to the microscope. High speed videos are captured 
using a Kodak Motion Corder camera and transferred to 
computer. 
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III. MYELIN GROWTH: THE FLUID FLOW OF 
SOAP 

Neglect for the moment, how a myelin forms in the 
first place. In this Section, let us consider what causes a 
myelin, once formed, to grow. There is no reason a priori 
to suppose the questions of how myelins form and how 
myelins grow can be separately considered; indeed, as I 
will suggest in Section III, myelin formation and myelin 
growth arise from the same phenomenon - the fluid flow 
of surfactant driven by the gradient in water concentra- 
tion. Nevertheless, it turns out to be both convenient 
and natural to consider these two questions separately. 

As previously reported, myelins grown in three exper- 
imental set-ups described in Section II follow different 
time dependences |12j . In contact experiments, it has 
been well-documented that the length L of the myelin 
bundles grows as L (x t^^"^ (Fig.^i). By contrast, myelins 
in drying drop experiments grow linearly in time: L oc t 
(Fig. [sj:). Myelin growth in immerse- and-puncture ex- 
periments follow yet a third time-dependence which can- 
not easily described by any simple function (Fig. ^p). 
A successful understanding of how myelins grow should 
describe all these different behaviors. 

From the drying drop experiment, we know myelins 
grow by the flow of surfactant from the root; from the 
contact and immersion experiments, we know applying 
osmotic stress can suppress myelin formation [T^]. For 
myelin-forming surfactants, the surfactant chemical po- 
tential /is decreases sharply with increasing <i>, the local 
water fraction in the lamellar structure, reaching a mini- 
mum when $ = <i>cq « 0.5, the fully swollen value [16 . If 
the surfactant sample is not uniformly hydrated, the hy- 
dration gradient V/is[<I>(x)] is a thermodynamic pressure 
difference; here /Xg is the surfactant chemical potential, 
and <i>(x) the local water fraction in the surfactant sam- 
ple. Provided the surfactant molecules are mobile, the 
hydration gradient can drive a surfactant flow from re- 
gions with low water content to regions of higher water 
content. In the simplest case, one can model the volu- 
metric surfactant flow as 

Js = -7sVms[*(x)] « 7s/i:'^V$(x) = aV$(x) , (1) 

where 7s is the surfactant mobility and fj,^^ ~ fee T is the 
fuUy-hydrated surfactant chemical potential; cr = 7s/is'^ 
is therefore a diffusion constant. Myelins result when this 
flow breaks through the initial surfactant/ water interface 
S, and the growth rate given by 

= Js • = crV$(x) • ; (2) 

here L is the myelin length, and n is the outward unit 
normal on S. Generally speaking, surfactant flow is not 
the only one present: water can flow too. The hydration 
profile ^(x) can evolve in time, for example, as water 
permeates into the dehydrated surfactant, 

^ = ^V^<l>(x,t); (3) 




t[s] 

FIG. 5: Myelin growth has different behaviors in different 
experimetal geometries, (a) Contact experiments: solid line 
is a power law fit: L oc t^^^. (b) Immersion-and-punction 
experiments: solid line is a fit to Eq. [9] insert shows the same 
data and fit plotted on logarithmic axes, (c) Drying drop 
experiments: solid line is a linear fit: L oat. 



here D is the diffusion constant for water permeation. 
The full dynamics of myelin growth arises from the cou- 
pling of Js to <i>. In particular, since transverse surfactant 
exchange between adjacent but unconnected bilayers is 
negligible in myelin-forming surfactants [17], the surfac- 
tant diffusion constant a — (c||,cr_L) is a tensor whose 
transverse component a± is zero; therefore surfactant 
flow is possible only if the hydration gradient has a lateral 
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component. As I will show below, once the proper bound- 
ary and initial conditions are imposed, this minimal set 
of ingredients are sufficient to reproduce quahtatively all 
the observed myelin growth behaviors. 

It is useful here to have an estimate of the values of the 
transport coefficients a\\ and D. The lateral diffusion of 
surfactant molecules on bilayers has been well-studied. In 
the limit of small driving forces, one can identify crj| with 
the bilayer's lateral surfactant diffusion constant: a\\ ~ 
2 — 3 /Ltm^/s. On the other hand, how water moves within 
a surfactant lamellar structure has not been studied in 
detail, and I rely instead on the following estimate: if the 
movement of water within a surfactant sample is rate- 
limited by the permeation of water across bilayers, then 
in the limit of low forcing, D ^ Kd^;, where K is the 
bilayer permeability and c?w is the inter-bilayer separation 
[IHllIilin]- For DMPC, if « 2 fim/s, ~ 3 nm, yields 
D « 6 X 10-3 /imVs. 



A. Contact Experiment 

In the contact experiment, myelins form at the ex- 
tended interface between a large mass of dry surfactant 
and bulk water (Fig.|2|. The situation may be simplified 
as a semi-infinite domain of surfactant (x < 0) contacting 
a semi-infinite domain of water (x > 0), with the inter- 
face at a; = 0. The problem reduces to a system of two 
ID equations. 



— =D— and — = a — 
dt dx'^ dt dx 



(4) 



x=0 



governing water equilibration and myelin growth, respec- 
tively. Here, I have replaced a with a coarse-grained av- 
erage (T, since in these experiments, the surfactant sample 
is always very disordered, with small randomly oriented 
lamellar domains [T2]. The initial and boundary condi- 
tions are 



B. Immersion Experiment 

In the immersion experiment, initially the hydration 
gradient is transverse to the bilayers of the well-ordered 
lamellar structure. However, the transverse surfactant 
diffusion constant is crj_ = 0; without surfactant flow, 
myelins do not grow |12j . When the needle is pushed 
into the lamellar structure, the puncture acts as a point 
conduit for water invasion into the surfactant (Fig. [sja) . 
There, the hydration gradient will be radial-lateral in ori- 
entation, to which the lateral surfactant diffusion (con- 
trolled by cry) can couple, leading to myelin growth. One 
can approximate this situation as a ID problem in polar 
coordinates. The relaxation of the hydration profile and 
myelin growth are given by 



(9$ 



= D 



1 a$ 

r dr 



Qj.2 



, dL 

^^^^ T^r = ^11 ^ 

dt " or 



, (7) 



with the initial and boundary conditions 



$(r > ro, 0) = $(c3o, t) = , $(ro, t) = $eq ; (8) 

here 2ro is the effective size of the water conduit created 
by the needle. Solving for L{t) yields j2Ij 



m = 



Jo [J^,{u)+Y^{u)] 



du 



(9) 



Compared to data, Eq. [9] gives a reasonable fit 
(Fig. [5|d). Again, the solution contains a prefactor pro- 
portional to <i>eq, indicating that myelin will be sup- 
pressed by applied osmotic stress. From the fit to Eq. [9] 
one obtains D/r^ k, 0.01 s^l and a\\rQ/D « 200 ^m. If 
we take cru = 2 /im^/s, then the fit implies ~ 1 /im as 
the effective size of the puncture, a small but physically 
plausible value. The fit also implies D « 10"^ ^m^/s 
as the diffusion constant for water permeation; this is 
consistent with the independent estimate of D based on 
bilayer permeability. 



^{x < 0, 0) = $(-cx), = , $(0, t) = $oq 
The solution for myelin growth L(t) is given by 

2a$e 



2q ^1/2 



(5) 



(6) 



which has the same time-dependence as observed in ex- 
periment (Fig.js^i). The prefactor in front of t^^'^ is pro- 
portional to the equilibrium water volume fraction <i>oq. 
Under the application of osmotic stress H, <I>eq <I>eq(H) 
is suppressed; this implies myelin growth will also be sup- 
pressed, in qualitative agreement with experiments |12j . 
Using CT « CT|| = 2 /im^/s, D = 6 x 10-^ ^m^/s, and 
<i>oq — 0.5 into the prefactor m, one obtains 15 /im-s"^/^ 
for the prefactor, which is consistent with experimentally 
measured values of 20 — 40 ^va ■ s~^/'^. 



C. Drying Drop Experiment 

In the drying drop experiment, myelins arc formed 
from pancakes situated at the pinned contact line of a 
slowly evaporating drop (Fig. |4|. Because the contact 
line is stationary, water influx across the contact line into 
the pancake must be balanced by water loss due to evap- 
oration. The hydration profile <&(x) is therefore static, 
and dL/dt cx V<i>|s is a constant. Myelin growth should 
be linear in time, 



L{t) = anV^ls t 



(10) 



in agreement with experiment (Fig. [5p). Fitting the ex- 
perimentally observed L(t) yields myelin growth rates of 
0.2 — 0.5 /im/s. This implies at the strength of the hy- 
dration gradient V<i>|5 is approximately 0.1 — 0.3 /im^^ 
at the contact line, which is physically plausible. 
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FIG. 6: (color online) (a) The deformation of a pancake to form a myelin. A is the pancake area (shaded); P, L are respectively 
the pancake perimeter and the myelin length (traced in thick outline); the myelin forms a,t t — to. (b) The width w at the 
waist of the elongated pancake as a function of to — t. 



IV. MYELIN INSTABILITY: FLOW AND 
BUCKLING 

The drying drop experiment provides the clearest pic- 
ture of how a myelin is formed: the initial structure is 
well-ordered, the dynamics are slow, and the view is un- 
obstructed. In this experiment, isolated, slowly growing 
myelins are formed from small, planar surfactant "pan- 
cakes" situated at the pinned contact line of a slowly 
evaporating water drop. The process follows a charac- 
teristic sequence: the initially rounded pancake becomes 
elongated, develops a pronounced "waist" which quickly 
narrows down to a cylinder — the myelin (Fig. [6^, in- 
sert). The entire event take place in ~ 10 s, at the end of 
which the original pancake has transformed into a pan- 
cake with a myelin growing from its leading edge. 

By tracing the pancake/myelin outline with a multi- 
point spline, the pancake morphology can be easily mon- 
itored as the myelin forms. Although the formation 
of a myelin represents a drastic change in morphology, 
the area A of the pancake remains constant through- 
out. Meanwhile the combined boundary length P + 2L, 
where P is the perimeter of the pancake alone and L 
the myelin's length, smoothly increases across the tran- 
sition (Fig. [6^). This suggests myelin formation may 
be thought as the deformation of a 2D shape (the pan- 
cake) of fixed area subjected to a continuously increasing 
boundary, with the myelin being made up of two lengths 
of excess perimeter. 

There is a surprising resemblance between myelin for- 
mation and the break-up of a pendant liquid drop. Both 
start with a rounded shape which elongates and develops 
a rapidly narrowing waist — although the drop eventu- 
ally breaks-up while the myelin stays attached. In viscous 
drop break-up, the radius r of the waist narrows in time 
as r oc to — t |22j . where to is the instant of break-up. In 
myelin formation, using high speed (250 images/s) video 
microscopy, the waist w(t) of the pancake is found to 



narrow in a similar way, 

w-doc (to-t)'-'*°-\ (11) 

(Fig. [6|d); here d is the myelin diameter, and to is the 
instant when the myelin formed. 

The resemblance between myelin formation and drop 
break-up, I suggest, is not coincidental but reflect similar 
underlying physics. The pendant drop is a 3D object of 
fixed volume. Gravity stretches the drop, surface tension 
resists the deformation, and the competing forces lead to 
a Rayleigh instability, breaking up the drop. The pan- 
cake is a quasi-2D object of fixed area. It is stretched by 
surfactant flow, bilayer elasticity resists this deformation, 
and the competing forces lead to mechanical instability, 
forming a myelin. As I hope to show, in myelin formation, 
as in drop break-up, the competition between a deform- 
ing body force and a restorative boundary force leads to 
an instability and rapid morphological change. 

Here, the mixed mechanical properties of the surfac- 
tant lamellar structure are crucial. A bilayer is a fluid 
membrane within which surfactant molecules can freely 
diffuse; but the bending the bilayer costs elastic energy 

^''^^JJ(Ci + C2fdS; (12) 

here k (^ 20 fceT") is the bending modulus, Ci, C2 the 
principle curvatures, and the integral is over the entire 
membrane. Stacked into a lamellar structure, the bilayers 
maintain a preferred repeat spacing, resisting compres- 
sion and dilation; moreover, for poorly-solubale myelin- 
forming surfactants, material exchange between discon- 
nected bilayers is negligible [T7|. Thus laterally, parallel 
to the bilayers, the lamellar structure acts as a fluid; but 
in the transverse direction it is an elastic solid. 

A pancake astride the pinned contact line of a water 
drop is more hydrated at its leading edge, inside the drop, 
than at its rear. As described in Section III, we can 
describe myelin growth on the basis of surfactant flow 
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FIG. 7: (color online) (a) The pancake/myelin junction, (b) Local balance of forces determines the static junction angle 6. (c) 
The distribution of 6 observed in experiment, (d) The combination P + L as a function of time; P is the pancake perimeter, L 
is the myelin length. P + L is clearly discontinuous at to, when the myelin forms. This is in contrast to Fig. |6^, which shows 
P + 21/ is continuous across to. 



driven by the hydration gradient across the pancake; the 
same flow of surfactant from the rear of the pancake to- 
wards its leading edge can also stretch and deform the 
pancake, leading to myelin formation. Additionally, such 
a flow also increases the overall hydration of the pancake. 
As Huang et al. showed, even a slight increase in overall 
hydration can overcome the added elastic energy of form- 
ing a myelin [23]. In this sense, the hydration gradient 
plays the role in myelin formation analogous to the force 
of gravity in pendant drop break-up. 

Confocal microscopy indicate the pancake is composed 
of multiple nested, closed bilayers enclosing little free vol- 
ume [12 . Therefore even a flat pancake has some bend- 
ing energy because the bilayers are sharply folded at the 
pancake's edge (Fig.|8|L). Since the bilayer curvatures are 
localized at the pancake's edge, and since the pancakes 
are thin (1 — 2 /im), a pancake's total bending energy 
is just proportional to its perimeter, — XpP; Ap is 
then the pancake's line tension, the 2D analog to the 
surface tension of a liquid drop. We can see its effect in 
the pancakes' rounded shapes prior to myelin formation. 
Presumably, Ap also acts to resist the deformation of the 
pancake caused by surfactant flow. 

To complete the analogy with pendant drop break-up, 
what instability, if any, does the competition between 
extensional surfactant flow and pancake line tension lead 
to? There is no Rayleigh instability in 2D fluid flow, 
however, the geometry of the myelin/pancake junction 
(Fig.|7|i) hints at the existence of a mechanical instability. 
The pancake has a line tension Ap; similarly, the myelin 
has a line tension (bending energy per unit length) Am. 
Local force balance at the junction between pancake and 
myelin leads to the relation 

Am = 2ApCos(0/2) , (13) 

where 9 is the static junction angle (Fig. [7|d). Eq. [T3]has 
a solution for 9 only if A„i < 2Ap, i.e. forming a myelin 
from two lengths of perimeter leads to a reduction in 



bending elastic energy. Experimentally, 9 is measured 
to be 110°-130° (Fig. which implies A^ ~ Ap. If 
this implication is true, then the total elastic energy of 
the pancake plus myelin is proportional to the pancake 
perimeter plus myelin length, 

.F^+m = ApP + Xn,L (X P + L , if Ap « Am . (14) 

As a function of time, P -I- L is discontinuous at to, the 
instant of myelin formation: it drops sharply, before re- 
suming growth at a slower rate (Fig. |7|i). This suggests 
^p+m is discontinuous at indicative of an instability. 

The instability identified here is related to the collapse 
of an elastic tube. The transverse cross section of an 
elongated pancake, as sketched in Fig. [8^, should be the 
same as that of a collapsed elastic tube {e.g. a soft wa- 
ter hose when the pressure head is turned off). This is 
because transverse cross sections of (thin-walled) elastic 
tubes are closed curves with minimal bending energies 
|24j . Starting with a unperturbed tube with a circular 
cross section, by lowering the pressure p inside the tube, 
the tube will eventually collapse: the cross section will 
have a flat region of wall-to-wall contact in the middle, 
and teardrop-shaped bulges at the edges [33 The 
shape of the transverse cross-sectional is controlled not 
by internal pressure per se but the ratio e = c^/a, where 
c is the cross-sectional circumference, and a the cross- 
sectional area. For an elastic tube, c is fixed; what low- 
ering p does is to reduce a. This is how the tube cross 
section evolves from a circle, where e is minimal, to a 
collapsed shape, where e is large. 

In myelin formation, the constraints are reversed. 
Here, because the pancake has a fixed thickness (lamellar 
structure is a solid in the transverse direction) , the area a 
enclosed by a transverse cross section must be fixed. In- 
stead, because a pancake is laterally fiuid and can flow, 
as it elongates, the width w — and circumference c — of 
a transverse cross-section in waist area will decrease. For 
a sufficiently elongated pancake, its width at the waist 
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FIG. 8: (color online) (a) Schematic drawing of a pancake/myelin complex, with transverse cross sections (shaded) showing 
its nested lamellar structure. The transverse cross section of the pancake looks like that of a collapsed elastic tube, while the 
transverse cross section of the myelin is circular, (b) The cross-sectional bending energy as a function w, the width of the 
transverse cross-section; d is the diameter of the resulting myelin. Shown above the curve is a sequence of calculated cross 
sections that would be obtained by a transverse cut across the narrow waist of the elongated pancake, as indicated by the 
arrows in the photographic insert. These show the transition from elongated pancake (wide cross section that is collapsed in 
the middle) to the myelin (circular cross section); the transverse cross section narrows in width approaching myelin formation, 
but the area a it encloses is constant. 



can become too narrow (e too small) to sustain a col- 
lapsed shape, and the circular transverse cross section of 
a myelin will become favored. 

In the linear elastic regime, the family of transverse 
cross sections from pancake to myelin (Fig. ^p), where c 
is variable but a is fixed, can be obtained numerically; 
I follow the calculation of Flaherty et al. for collaps- 
ing elastic tubes [25], modified to reflect the constraints 
appropriate to myelin formation (see Appendix). The 
bending energy of each transverse cross section 



one can construct only a single dimensionless group. 



C'ds 



(15) 



(here C is the curvature, and the integral is along the 
cross-sectional circumference) can then be evaluated as a 
function of the cross-sectional width w. Plotted against 
ui, flattens to a constant (= 2Ap) when w is large; but 
when w < 2.3(i, the collapsed cross section is unstable: 
falls with decreasing w, reaching its minimum value 
(= Am) s.t w = d and a circular cross section — a myelin 
(Fig.lSlD). This calculation indicates Am « Ap, consistent 
with the estimate made on the basis of balance of forces 
at the myelin/pancake junction. Inserting the calculated 
value of Am/Ap into Eq. 13 yields a myelin/pancake junc- 
tion angle = 119°, in good agreement with observation. 

If the instability described above is indeed responsible 
for myelin formation, then its dynamics should be con- 
trolled by surfactant flow and bilayer elasticity. Just as 
dimensional analysis can be used to analyze the dynamics 
of fluid instabilities, e.g. drop break-up, a similar analysis 
can be applied to myelin formation. The relevant physi- 
cal parameters here are the 2D pancake viscosity -q-p (di- 
mension MT~^), the line tension Ap (dimension MLT~^), 
the reduced waist width w — d, and Iq — t. From these 



A 



p 



tp-t 
w — d 



d (xto — t 



(16) 



The implied linear dependence of w on tQ — t is consis- 
tent with experimental observations. One may further 
estimate the velocity Ap/ryp. Taking Ap — Am ^ Nn/d, 
where N w 100 is the number of nested bilayers, k w 
10^^^ J is the bending modulus of a single bilayer, and 
d « 1 /im is the myelin diameter; and taking -q-p ~ Nrj, 
where 77 « 10 cP • /xm is the 2D viscosity of a single bilayer 
[28], one obtains Xp/rjp ~ 10 /im/s. This compares well 
to the experimentally measured velocity (w 1 fim/s). 

Finally, this model suggests the hydration gradient is 
not the only force which can lead to myelins. Any body 
force (such as shear stress) can also lead to myelin for- 
mation, as long as it can sufficiently stretch the pancake. 
Indeed, myelins can be induced to form by shear flows 
even in the absence of hydration gradients [T2] . 



V. SUMMARY 

In this paper, I argue the fluid flow of surfactant, driven 
by the hydration gradient, results in the formation and 
growth of myelins. The immediate consequence of this 
argument is there can be no myelins without surfactant 
flow. This give a natural explanation to the well-known 
observation that myelins do not form for temperatures 
below the so-called main transition temperature Tm (0°C 
for DLPC and 23°C for DMPC) [22 . Below Tm, the 
surfactant tails freeze into a solid phase and the lamellar 
structure loses its lateral fluidity; consequently myelins 
can no longer form. 
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The varying growth behaviors as seen in the different 
experiments, I suggest, is entirely due to the different 
ways the driving force, the hydration gradient, relaxes 
in a given situation. This is a function of geometry, as 
well as initial and boundary conditions. A good test of 
the model suggested in Section III is to devise additional 
experimental geometries, and to examine myelin growth 
in those instances. 

While the fluid flow of surfactant is crucial, I suggest 
the initial formation of myelin figures also depends on 
the elastic character of the surfactant lamellar structure. 
In analogy to a classic fluid instability — the break-up 
of a pendant drop — myelin formations is a result of the 
competition between bulk surfactant flow, driven by the 
hydration gradient, and bilayer bending, manifesting it- 
self as a line tension. The instability at its heart is closely 
related to a classic elastic instability — the buckling and 
collapse of an elastic tube (in reverse) . The elastic forces 
involved here may be measured directly. By pulling on a 
well-ordered, equilibrium surfactant lamellar pancake un- 
til a myelin is formed, using possibly an optical tweezer 
or a microcantilever, one should be able to determine the 
pancake edge tension Ap and the myelin line tension Am 
from the work of deformation. This would be a direct 
test of the myelin instability mechanism suggested here. 

The experiments and calculations described here aim 
to show that the formation of myelin depends crucially on 
the mixed mechanical properties of the surfactant lamel- 
lar structure: the wet soap has to be able to flow as a fluid 
and buckle like an elastic solid. Mechanical instabilities 
in fluids {e.g. the Rayleigh instability) and elastic solids 
{e.g. buckling) have been active subjects of study since 
the work of Euler; myelin formation is of a rather less fa- 
miliar type: a mechanical instability of mixed character 
in a material where fluid and elastic properties coexist. A 
large number of similarly mixed instabilities must surely 
exist, given the diversity of liquid crystalline phases, both 
thermotropic (temperature controlled) and lyotropic (as 
in soap). In particular, some of the diverse interfacial 
structures seen in surfactant dissolution may be under- 
stood in terms of mechanical instabilities of the interface. 

Finally, I have not treated an especially striking in- 
stability of the myelins themselves: coiling and the for- 
mation of helices (see Fig|l]) . The coiling often appears 
spontaneously, and can also be induced by the addition 
of polymers or Ca^+ ions [S [M 123 ED]- The case of 
polymer-induced coiling has been attributed to the ap- 
pearance of a spontaneous membrane curvature in the 
presence of the additive. But in general one expects that 
myelins, being laterally fluid, cannot support the tor- 
sion needed to form helices and coils. The treatment of 
myelins in this paper is simplistic in the sense I assume 
the myelin is perfectly fluid laterally and perfectly elastic 
in the transverse direction. While such a simple picture 
appears sufficient when it comes myelin formation and 
growth, the coiling instabilities requires more elaborated 
treatments [5T1 [5^ . 



APPENDIX A: CALCULATING THE ELASTIC 
ENERGIES 

1. Linear elasticity: thin-walled elastic tube 

The calculation for the transverse cross section of pan- 
cakes and myelins is based on the calculation of Flaherty 
et al. [2S] for collapsed elastic tubes; it thus useful to 
briefly outline the physics which gives us the cross sec- 
tions of elastic tubes. 

The calculation of Flaherty et al. assumes linear elas- 
ticity. For a tube made from an conventional elastic ma- 
terial, linear elasticity holds if strains are small; practi- 
cally, this is satisfled if the tube wall is thin. Then, the 
equations of mechanical equilibrium, which one solve to 
flnd the tube cross sections Xt(s), can be derived from 
the theory of thin elastic shells ^B]- They can also be 
derived from minimizing the energy functional 

7^[xt(s)] = | ^(C[xt]-l/rt)'ds-| jn-^tds- (Al) 

here C is the radius of curvature, rt is the radius of the 
unperturbed tube, p is the pressure difference between 
the interior and exterior of the tube, and fi in the outward 
unit normal [53]. One seek solutions xt(s) of the form 



xt(s) 



f 


cos6(s) 


'0 


sin 0(s) 



ds 



C{s) ^ e'{s) ; (A2) 



here c = 27rrt is the tube's cross-sectional circumference. 
Since c is flxed for an elastic tube, one further requires 



xt(s) = xt(s-(-c), e(s + c) = 6 + 27r 



(A3) 



The second term in 7i, involving p, is simply —pa, 
where a is the area enclosed by xt. This term constrains 
the tube cross section, whose circumference c is fixed, to 
enclose an area a (which is set by the pressure difference 
p). The first term in Ti, involving the curvature C, is the 
elastic energy of the transverse cross section; this term is 
just f^, the cross-sectional bending energy as defined in 



Eq. 15 plus a constant. Thus xt is a cross section of cir- 
cumference c (which is fixed), enclosing an area a (which 
varies as a function of p) , whose cross-sectional bending 
energy is minimal. Finally p, being a dimensional 
quantity, cannot be the parameter which characterizes 
the cross sectional shape; instead, that role is played by 
the dimensionless ratio e = c^/a. 

The transverse cross sections of the pancake, like those 
of the elastic tube, minimizes the cross-sectional bending 
energy for a given e. But here the constraint to shape 
evolution is different: a is fixed and c varies as cross- 
sectional width w varies. The control parameter is still 
£, and so we can obtain a pancake cross section simply by 
re-scaling a similarly shaped elastic tube cross section. If 
Xt (s) is an elastic tube cross section with fixed circumfer- 
ence c enclosing a variable area a, a pancake cross section 
Xp(s') of the same shape (having the same value of e). 
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with variable circumference d enclosing a fixed area a', 
is given by 



xp(s') 



1/2 



xt(s) . 



(A4) 



Once family of elastic tube cross sections has been found, 
applying Eq. |A4| will yield a corresponding family of 
pancake cross section, whose widths and cross-sectional 
bending energies can be evaluated. The curve in Fig. 4b 
was calculated from 1500 pancake cross sections, whose 
the cross sectional width w was varied systematically 
from Ix to ^ 4x the myelin diameter d. 



2. Thick walls: effective bending modulus Keff 




FIG. 9: (color online) The transverse cross section of a pan- 
cake with wall-thickness I. The solid line is the median surface 
s; the dashed line s(/i) is a surface (bilayer) displaced by h 
from s. 



The agreement between the previous calculation, based 
on thin-walled tubes, and experimental observation of 
myelin formation may be surprising. With a thin-walled 
tube, any radius of curvature will be large compared to 
wall thickness, therefore strains ares small and linear elas- 
ticity applies. With thick-walled tubes there can be large 
strains and nonlinearities. 

For a conventional elastic material, if a slab of thick- 
ness i is bent through a radius of curvature i?, the max- 
imum strain within the slab will be on the order of £/-R; 
when i? ~ £, the maximum strain is of order 1 and the lin- 
ear elasticity breaks down. This means the slab's bending 
modulus will be function of the radius of curvature rather 
than a constant. 

It may appear that the previous calculation cannot ap- 
ply to a thick-walled pancake with N ^ 100 bilayers. 
However, a lamellar structure of surfactant bilayers is not 
a conventional elastic material: the bilayers are fluid, and 
adjacent bilayers are free to slide past each other. Thus 
the bending of each bilayer can be considered indepen- 
dently, subject to only the requirement of a fixed bilayer 
repeat spacing do; the total cross-sectional bending en- 
ergy of a iV-layer pancake is just the sum of each bi- 
layer's individual contribution. We can rewrite this sum 
in terms of the bending of a single median surface s and 



a curvature-dependent effective bending modulus Kcfj 

N 

= f E / ^'i^'- ^If «cff^2ds . (A5) 

For large N, we can make the continuum approximation 
where the bilayer s,; — > s{h) is labelled by its displace- 
ment h from the median surface. 



e/2 



ds{h) dh 
-e/2J {R+hfdo 



Kcff 

i?2 



ds 



(A6) 



here £ = A^do and R—1 jC is the radius of curvature at 
the median surface (Fig. [9]). 

If the curvature C on s is approximately constant 
over length scales comparable to £, then due to the fact 
the lamellar repeat spacing is fixed, the differential arc 
lengths ds{h) on surface s{h) and ds on s are related by 



ds{h) 



R+h 
R 



ds 



(A7) 



Inserting this relation into Eq. |A6| yields the following 
expression for KcS, 



i?\ /R + £/2 

KeS = K \ — ] log 



da 



R-e/2 



(A8) 



As shown in Fig. |10[ HqS has significant variations with 
R only close to the singular limit oi R = £/2. Therefore 
a constant bending modulus independent of local curva- 
ture, which is a consequence of linear elasticity, should re- 
main an excellent approximation. The analysis for myelin 
formation based on thin-walled tubes should hold even 
for thick- walled pancakes with many bilayers. 



^ 2 

1 



R/i 



FIG. 10: The effective bending modulus Keff of A'^ bilayers as 
a function of R, the median radius of curvature; £ — Ndo is 
the total thickness of the A^'-bilayer lamellar structure. 
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